Quantum Hall effect in a singly and doubly connected 3D topological insulator 
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The surface states of topological insulators, which behave as charged massless Dirac fermions, are 
studied in the presence of a quantizing uniform magnetic field. Using the method of D.H. Lee^, 
analytical formula satisfied by the energy spectrum is found for a singly and doubly connected 
geometry. This is in turn used to argue that the way to measure the quantized Hall conductivity 
is to perform the Laughlin's flux ramping experiment and measure the charge transferred from the 
inner to the outer surface, analogous to the experiment in Ref.— . Unlike the Hall bar setup used 
currently, this has the advantage of being free of the contamination from the delocalized continuum 
of the surface edge states. In the presence of the Zeeman coupling, and/or interaction driven 
Quantum Hall ferromagnetism, which translate into the Dirac mass term, the quantized charge Hall 
conductivity a xy = ne 2 /h, with n = 0,±1,±3,±5. . .. Backgating of one of the surfaces leads to 
additional Landau level splitting and in this case n can be any integer. 
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Theoretical prediction of the existence of an odd num- 
ber of Dirac cones in the dispersion of the surface states 
of topological insulators^ and subsequent experimen- 
tal observation of such unusual surface states^— has pro- 
pelled this research field into an active area of condensed 
matter physics (for reviews see Refs^— ). Particularly 
interesting is the problem of the topological insulator 
surface Dirac Fermions in magnetic field. Because the 
Dirac Fermions carry definite charge, the magnetic field 
couples to the orbital motion. If this motion is con- 
strained to be perpendicular to the applied field, the Lan- 
dau level quantization results. However, as discussed by 
D.H. Lee 1 , since the Dirac Fermions move on the sur- 
face of a 3-dimensional material, in the absence of mag- 
netic monopoles, i.e for V • B = 0, it is impossible for 
the magnetic field to be everywhere along the normal 
of an oriented surface. Instead, in a typical experimen- 
tal setting, the three dimensional material is placed in 
a uniform external magnetic field, and only portions of 
the surface, say the top and the bottom ones, experience 
Landau quantization. The Dirac Fermions on the sur- 
faces tangential to the external magnetic field continue 
moving as plane- waves. 

In addition, the spin-orbit coupling, which causes the 
appearance of the Dirac particles in the first place, makes 
the Zeeman coupling different from that in graphene, 
where Dirac particles also appear but where the spin- 
orbit coupling is negligible. Thus, instead of simply spin 
splitting the electronic energy levels, the Zeeman term 
in topological insulators acts as a Dirac mass. As illus- 
trated in FigUH this causes the splitting of the zeroth 
Landau level, but the higher Landau levels are not split 
unless their guiding center approaches the edge. Rather, 
at positive energies they move up and at negative energies 
they move down. Of course, because of the Dirac struc- 
ture, the energy scale associated with Zeeman splitting 
~ 1 K x H[T] is much smaller— that the spacing between 
the zeroth and the first Landau levels, ~ 200^ x y/H[T] 
for realistic fields. Nevertheless, similar behavior can be 
expected in the presence of interaction driven Quantum 



Hall ferromagnetism where the "Zeeman" scale can be 
much larger. While neglecting orbital coupling, the ef- 
fects of Zeeman splitting near the boundary between the 
portions of the material where the field is perpendicular, 
and therefore the Dirac point is gapped, and where it is 
parallel, and therefore gapless, was also studied in RefJ^. 
Some aspects of Landau quantization in thin films of 
topological insulators were also analyzed in Ref J^. The 
analytical results for the energy spectrum obtained in 
this work are in agreement with recent numerical lattice 
band-structure diagonalization studies in magnetic field 
for systems with small withdthi^. 

For a sphere with a finite radius r the Dirac Hamil- 
tonian adopted to this curved surface in the presence of 
a uniform applied field (without Zeeman coupling), was 
analyzed in an insightful study in Ref*i. In the Lan- 
dau gauge, the eigenenergy is a function of the azimuthal 
quantum number m which serves as a "guiding center" 
coordinate. For large and positive values of m the spec- 
trum corresponds to nearly doubly degenerate Landau 
levels, with the wavefunctions residing near the north and 
south poles. The energy splitting is exponentially small 
in ~ t/Ib% where the magnetic length £b = \Jh~c/eB. 
As m approaches zero, the wavefunctions move towards 
the equator of the sphere, the Landau levels are split 
and merge into the plane-wave-like states residing near 
the equator. The spectrum of these states resembles (fi- 
nite size) Dirac spectrum. Thus, for the singly connected 
topology, the portion of the surface approximately tan- 
gential to the external magnetic field harbors the chiral 
quantum Hall edge states together with the conducting 
non-chiral surface states, which disperse as Dirac parti- 
cles. Since the latter do not localizei^, the two terminal 
conductance is not quantized even when the Fermi en- 
ergy lies between Landau levels. It was proposed^ that 
the way to measure the quantized Hall effect is to set up 
a potential difference between the electrode placed in the 
caps of the sphere and to measure the circulating current 
on the outer surface. 

Experimentally, quantum oscillations originating from 



the surface states have been reported by several 
groups 1 ^—. In each case (although to varying degree) 
the signal is contaminated due to the finite 3D bulk con- 
ductivity. In 70nm thick strained films of HgTe, quantum 
Hall effect has been reported 21 , although interestingly 
the quantized Hall plateaus appear without the longitudi- 
nal resistance R xx reaching zero. So far, no quantization 
plateaus of the 2D Hall conductivity has been reported 
in thick samples. It is important to ask whether such 
quantization could be realized in a typical multi-terminal 
experimental setup even for a perfectly insulating bulk. 
Theoretically, it was argued in RefJ^ that for topologi- 
cal insulators with odd number of Dirac points, the four 
terminal Hall conductance should remain quantized even 
in the presence of scalar disorder, although the six ter- 
minal should not. However, even in the ballistic limit, 
such quantization of the four terminal Hall conductance 
depends sensitively on the number of non-chiral surface 
channels, which may vary between the four electrodes. 
As such, it is sensitive to the surface roughness and for 
typical Fermi moment a 19 ' 21 kp ~ 0.2 — 0.5nm _1 would 
require few nm precision in the height of the sample, as 
in the molecular beam epitaxy grown films of Ref. 21 . 

Here we study both the singly and the doubly con- 
nected geometry. In the former case, the standard Hall 
bar setup (FigEJ) will not lead to quantization of the Hall 
conductivity, unless the sample height is reduced to be 
much smaller than the magnetic length. The latter case is 
illustrated in Fig[JJand argued to be an alternative way to 
measure the quantization of a xy . The idea is to perform 
the analog of the Laughlin thought experiment, experi- 
mentally realized in Ref^, and to measure the amount 
of charge AQ transferred from the inner surface to the 
outer surface in response to the change in the magnetic 
flux A<& threading the sample. Then, 

This setup has the virtue that any interaction driven 
fractional quantum Hall effects can also be detected, as 
shown in the context of 2DEG heterostructures in Ref. 2 . 

We build on the formalism and findings presented in 
Refi and analytically study the energy spectrum such 
geometry. We include the Zeeman coupling, as well as 
a difference in the gate voltage applied between the top 
and the bottom surface of the 3D sample. The former 
causes the splitting of the zeroth Landau level, but not 
the rest of the Landau levels, while the latter splits all 
Landau levels. Thus, integer quantum Hall conductivity, 
measured in the doubly-connected setup, can take on any 



positive or negative integer, including zero. 

This paper is organized as follows, in Sec. I the effective 
Hamiltonian for the geometry shown in FigJT] is derived 
and the boundary conditions at the edges are discussed. 
In Sec II, the matching of the wavefunctions at the cor- 
ners leads to the equation which determines the condition 
for the eigenenergies as a function of the guiding center 
k£s- This equation is shown analytically to describe a 
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FIG. 1: The physical setup proposed for the measurement of 

Dirac-like continuum as well as Landau levels. Its nu- 
merical solution determines the behavior of the discrete 
levels as the guiding center approaches the continuum. 
In Sec III the Hall bar and the Corbino-like geometry 
are further compared. 



I. HAMILTONIAN, EIGENSTATES AND THE 
MATCHING CONDITIONS 

The specific geometry considered is shown in FigJT] 
The magnetic field is assumed to be perpendicular to the 
surface of the hollow cylinder of inner radius R, shell 
thickness 26 and height 2a. As discussed by D.H. Leei, 
the Dirac Hamiltonian needs to be written on the sur- 
face of this two dimensional curved space. In the limit 
of R — > oo the system is equivalent to an infinitely long 
slab with rectangular cross-section and periodic bound- 
ary conditions along its axis. We can now use the polar 
coordinates shown in FigUJto describe the surface. 

The procedure for determining the Hamiltonian, which 
follows from the discussion in Refs.— >22, is shown the Ap- 
pendix. The result is 
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where periodic boundary conditions for x, which is al- 
lowed to get arbitrarily large, are assumed. Obviously, 
the wavefunctions separate and correspond to planewaves 
with wavector k in the ir-coordinate. This will serve as 
the analog of the guiding center coordinate. The term 
proportional to h z represents the Zeeman coupling. A 
mean-field description of any possible quantum Hall fcr- 
romagnetism would have a different value of h z in the 
first and the third region than in the second and fourth. 
As usual, the magnetic length is Ib — y/hc/eB. 

We can now find the eigenfunctions at a fixed energy 
E and a fixed k for each segment. This is done in detail 
in the Appendix. For the 1 st and 3 rd segment the Dirac 
eigenspinors can be written in a closed form in terms of 
the parabolic cylinder functions, D u (z), with indices v 
of the upper and the lower components differing by 1. 
For the 2 nd and 4" 1 segment, the eigenspinors are plane- 
waves. In what follows, we use dimensionless lengths and 
energy scales defined as 
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The matching of the wavefunction is discussed in de- 
tail below. While in general 4-parameter family of self- 
adjoined extensions is allowed, it is argued that for sharp 
edges, continuity of the wavefunctions should be im- 
posed. 



A. Matching conditions 

In order to completely specify the behavior of the wave- 
functions, the Hamiltonian ([2]) must be supplemented by 
boundary conditions at the points where the horizontal 
and the vertical surfaces meet. We will assume that the 
corners are sharp, meaning that the lengthscale associ- 
ated with the corner curvature is much smaller than the 
magnetic length £b- 

If we require that the Hamiltonian is a self-adjoint op- 
erator we require that any two spinor wavefunctions ipi 



and tp2 satisfy 

^\ L {M^2^2L{<t>Q) = ^\ R {M^2R{M, (5) 

where the subscripts L and R refer to the direction of 
approach of the boundary, i.e. left or rig h1j22. The most 
general linear homogeneous boundary condition imposed 
on ip2 is 



V>2L = Mtp 2 R, 



(G) 



where M is a 2 by 2 matrix. In order to satisfy ([S]) we 
then must have 

lpi L ((l>o)v2Mlp2R(<t)o) = ■0l fl (0o)CT 2 '02fl(0o)- (7) 

This must hold for arbitrary tp2R and therefore the 
boundary condition on tpi is 



(8) 



Since we require that the domains of H and H T coincide, 
the above must also be the boundary conditions on -02- 
Therefore, M must satisfy 

M- 1 = (j 2 M^a 2 . (9) 

Taking the determinant of both sides gives 



1 



= detM* => detM = e lx . 



detM 

Using the above relations we finally find that 
,. , W + Z X-Y 



X + Y W - Z 
1 = W 2 - Z 2 -X 2 +Y 2 . 



(10) 

(11) 

(12) 



That means that the self-adjoint constraint leaves us 
with 4 independent real paraments which determine the 
boundary conditions at the corners. The von Neumann- 
Weyl deficiency indices^ are therefore (2, 2) 

To determine these conditions we regularize the corners 
as small circles with the vertical and horizontal surface 
lines being tangential to the circle. Eventually, we are 
interested in taking the limit of the radius of the circle, 
ro, to zero. Note that in the limit tq <§; £b the vector 



II. THE ENERGY SPECTRUM 




FIG. 2: Energy spectrum for hz — and a — > oo. In addi- 
tion p — b/ls 3> 1. The shaded cones represent the Dirac 
continuum physically located at the inner (left) and the outer 
(right) surfaces of the doubly connected sample. Each Lan- 
dau level is doubly degenerate and at energy \/2nhvF / (-B- 
To understand how the discrete Landau levels merge into the 
continuum, see Fig lHIl 



potential does not depend on the position along the circle. 
As a result, the solution to the Dirac equation along the 
circle is a plane wave. The ratio of the two components of 
the Dirac spinor of the planewave solution is independent 
of the position, meaning that X = Y = Z = 0. At finite 
energies the wave vector must be finite and therefore in 
the limit of r —> the phase does not advance, i.e. % = 
0. This argument leads us to the requirement that for 
sharp corners, the Dirac spinors are continuous, or 

M = 1 2 . (13) 



In addition, we are going to take the limit b 3> £b and 
focus on the solutions with the "guiding center" k = kis 
either «6/^ = /3or« — /3. In the former case, K ~ (3, we 
can require that the wavefunction vanishes for tan <f> — > 
— oo on the top surface and for tan <f> — > oo on the bottom 
surface. Similarly, for k ~ — /?, we can require that the 
wavefunction vanishes for tan <j> —¥ oo on the top surface 
and for tan</> — ¥ —oo on the bottom surface. 



The energy spectrum at each k is determined from 
the matching conditions on the wavefunctions discussed 
above. Near each edge, these take the form of four lin- 
ear equations in four unknowns (|A37[) and (|A46p . with 
energy and k dependent coefficients. The technical de- 
tails are presented in the Appendix. Here we present the 
analysis of the energy spectrum which results from these 
equations. 

The non-trivial solution to the Eq. (IA37p exists pro- 
vided that 

[(w_ + i9-)B+ - eA+] [(w_ - i6_)A g+ - eB g+ ] e 2l9 - a = 
[(w_ - i9-)B+ - eA+] [(w_ + i6_)A g+ - eB g+ ] e - 2l9 - a . 

(14) 

The functions used in the above equation are derived 
in the Appendix. The above equation determines the 
energy spectrum for k ss j3, i.e. near the outer surface. 
If we are interested in the finite size effects, we have to 
keep a = ci/£b finite. This is useful to understand the 
degeneracy of each energy level caused by the tunneling 
across the vertical edge as well as mesoscopic transport 
effects discussed in the next section. On the other hand, 
in the thermodynamic limit a — > oo. 

In order to solve Eq. (fT4"l) for e = E j (yp j Ib) we need 
to consider two cases. 

For #2 > the left hand side of Eq.(fT4"f is the com- 
plex conjugate of the right hand side. Therefore, we are 
looking for roots of a purely imaginary function. Letting 

S = [(u- + i9_)B+ - eA+] [(u- - i9-)A g+ - eB g+ ] (15) 

we find that (1141) requires 

Note that the right hand side of the above equation does 
not depend on a while the left hand side has a single 
pole every time 29-a = mn where m = 0,±1,±2.... 
For a — > oo the spacing between the poles approaches 
zero and the left hand side changes from — oo to oo for 
mw/(2a) < 8- < (rn + l)7r/(2a). Since the right hand 
side is a function of e which changes on scales much 
larger than l/aasa->co, we have at least one positive 
and one negative eigenenergy solution for each interval 
m7r/(2a) < < (to + l)7r/(2a). This proves that for 

the energy spectrum forms a continuum as a — > oo. 

For 9 2 _ < 0, the left hand side vanishes in the limit 
a — > oo and the spectrum is determined by 

[(«_ + \9.\)B + - eA + ] [(w_ - \9-\)A g+ - eB g+ ] = 0. 

(17) 



For k — f3 <C —1, we can solve these equations within 
exponential accuracy by noting that for non-negative in- 
teger, the parabolic cylinder functions entering A+ and 
B + satisfy 



E 



D v (z) = 2-»/ 2 e-^H„(^=j , for i/ = 0,1,2,.. 



(18) 



where H n {z) is Hermite polynomial. If the v in D v (z) 
deviates from non-negative integer, the function diverges 
exponentially at large negative z. Therefore, as long as 
k — P <C —1, the equation 

[{LO- + \e_\)B + -eA + ] =0 

is solved for 

e = ±^2n + i] 2 z , n = 1,2,3,..., 

e = (19) 

Similarly, as long as k — (3 -C — 1 the solutions to 

[(u;_-|^_|M 5+ -eB 5+ ]=0 



arc 



n = 1,2,3,.. 



Vz 



(20) 



For k fa /3 the equations can be easily solved numerically 
and the dispersion of the Landau levels in the vicinity 
of the Dirac continuum can be determined. The solu- 
tion for -q z = v g = is shown in Fig|TT] for both edges. 
In this case all Landau levels are doubly degenerate and 
the a xy sequence is %- x (odd integer). Also, note the 
downward dispersion of the Landau levels as they ap- 
proach the continuum. For r\ z = 0.3, v g = the energy 
spectrum is shown in Fig. [TT| for outer the edge. The 
lowest Landau level is now split linearly in h z giving rise 
to the possibility that if the Fermi energy lies in this gap 
o xy = 0. The higher Landau levels are doubly degen- 
erate, but split as they approach the edge continuum. 
Therefore for v g = 0, the allowed values for a xy = ne 2 /h, 
are n — 0, ±1, ±3, ±5 . . .. Finally, for v g ^ all possible 
Landau level (double) degeneracies are lifted, and as a 

result n = 0,±1,±2,±3. . .. 

Similarly, the non-trivial solution to the Eq. (|A46[) ex- 
ists provided that 



[(u+ + i0+)B- - eA-} [(u+ - i9+)A g - - eB g -\ e 
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FIG. 3: Energy spectrum for hz ^ (r\ z = 0.3) and a — > oo 
near the outer surface. In addition (3 = b/ls 3> 1. The 
shaded cones represent the Dirac continuum physically lo- 
cated at the outer surface. The degeneracy of the Landau 
levels is split as k = kin approaches the continuum. The 
double degeneracy of the zeroth Landau level is lifted, with a 
gap that varies linearly in hz- The double degeneracy of the 
higher Landau levels is split near the continuum, but in the 
thermodynamic limit remains intact for k <<iC Mb with the en- 
ergies ±{Tivf /£) \/2n + 77I, n = 0, 1, 2, . . .. The spectrum near 
the inner surface can be obtained simply by taking k — > —n, 
i.e. by mirror reflecting the above picture around kls = 0. 
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The above equation determines the energy spectrum for 
k m — /3, i.e. near the inner surface. The analysis of this 
equation proceeds in the same way as for Eg. (1141) . The 
continuum appears for 

e 2 > ( K + f3 + r ]z ) 2 . 



(21) 



Furthermore, the equations for e 2 < (k + j3 + rj z ) 2 can 
be shown to be the same as (|14[) under the transforma- 
tion k — ¥ The spectrum near the inner surface can 
therefore be determined simply by mirror reflecting the 
spectrum near outer surface about k = 0. 



_kt B - b/t B 

-2.5 -2 -1.5 -1 -0.5 0.5 1 

FIG. 4: The positive energy spectrum for hz = and a — 
cl/Ib = 5 near the outer surface. In addition /3 = bjlB 3> 1. 
Note the presence of the discrete non-chiral modes, which 
merge into continuum for a — > oo. The spacing between the 
levels at his = fe/^s scales as ~ 1/a while the spacing be- 
tween the (Dirac) Landau levels for klB <C b/(.B scales as 
~ 1/£b- Therefore, in order to observe quantized Hall con- 
ductivity in the Hall bar geometry the condition a <C £b 
should be satisfied. 



III. DISCUSSION OF HALL BAR VS. 
"CORBINO" GEOMETRY 



Currently, the experimental geometry used to measure 
the Hall conductivity in the quantum Hall regime of 3D 
topological insulators is the Hall bar geometry sketched 
in Fig(5] No plateau quantization of a xy has been ob- 
served. While part of the reason for this is finite 3D bulk 
conductivity, we wish to argue here that even if the sys- 
tem was insulating in the bulk the presence of the non- 
chiral surface modes will spoil the quantization of a xy . 
One way to avoid such contamination would be to reduce 
the sample height <C Ib- The second way, would be to 
use "Corbino" geometry shown in Fig[TJ to ramp up the 
flux through the hollow region and to measure the charge 
transferred between the inner and the out surfaces. 

For the Hall bar geometry in the ballistic limit, the Hall 
conductance, as well as the longitudinal conductance, are 
easily obtained within the Landauer-Buttiker formalism. 
Assuming that the contacts 2, 3, 5 and 6 float to the av- 
erage chemical potential of the modes which enter them, 
for M-chiral modes and iV-non-chiral modes with perfect 
transmission the Hall conductance, measured between 



FIG. 5: (Top panel) A schematic for a Hall bar setup for the 
3D topological insulator. The applied magnetic field is per- 
pendicular to the top surface. The applied magnetic field is 
perpendicular to the top surface; six contacts are also marked. 
(Bottom panel) The corresponding edge/surface state struc- 
ture. There are M chiral modes coming from the top and 
the bottom surface Hall droplet. In the absence of the bot- 
tom gate M can be any odd integer (including zero if Zeeman 
field is included). In addition, there are N non-chiral modes 
coming from the surfaces parallel to the external field. 

contacts 2 and 6, is 

e 2 / iV 2 \ 
G Hall = j(M + N + -y (22) 

The longitudinal conductance measured between 2 and 3 
is 

e 2 / M 2 M 2 \ 

G = — [2N + M+ —— + — . (23) 

h V N N + MJ v ' 

For finite N, the Hall conductance here is generally not 
quantized. In the thermodynamic limit, N M and the 
Hall resistance is small (~ -j^M/N 2 ). These results were 
also obtained in RefJ^. 

The physical reason for the non-quantization for large 
sample height is quite clear. The presence of the large 
number of non-chiral surface states tends to make the lo- 
cal potential on each of the surfaces close to (/zl + Hr) /2, 
as opposed to [il on one, and [ir on the other, as is 
the case in the absence of non-chiral modes. In the 
presence of disorder, in the form of scalar potential 
in the Dirac equation, but in the absence of electron- 
electron interactions, the states have been argued to re- 
main delocalize d 15 ' 24 and the conductivity diverges as 
temperature T 0. Therefore in this limit there will 
be no voltage drop along the surfaces, all of which will 
appear near the contacts. The local potential near the 
two surfaces will still tend to (jil + [J-r.)/2 and no quan- 
tization of Gho.ii will occur. Including electron-electron 
interactions and scalar disorder, it has been argued in 
Refj2£ that the states remain delocalized, but that at 
T = the conductivity flows to a finite value. Therefore, 
the potential will drop along the surfaces, the modes will 



equilibrate but no quantization will occur. The two ter- 
minal conductance will also not be quantized^. The way 
to achieve quantization in the Hall bar geometry is there- 
fore to eliminate the side surface modes altogether, which 
can be achieved through finite size effects by reducing the 
height of the sample below £g. 

A four terminal setup shown in Fig. 6a of Ref.— has 
been argued to lead to quantized Hall conductance 
Me 2 /h even in the presence of disorder. However, this 
result was obtained assuming that the number of non- 
chiral channels is exactly TV for each of the four seg- 
ments of the four terminal setup. Unlike the number 
of chiral channels M, for a thick sample the number 
of non-chiral channels can vary from segment to seg- 
ment. In the four terminal setup, we should therefore 
consider Ni channels between \xl and p 2 , N 2 between 
p, 2 and 11 r, N3 between p,R and p^ and N4 between /X4 
and hl- In the ideal case of perfect transmission, we 
find h-h = T ( 2M + N 1 -N 2 +N 3 - N 4 ) (p L - fi 2 ). 
Therefore, unless Ni — N 2 = N4 — N3, the Hall conduc- 
tance defined this way is not properly given by the num- 
ber of chiral channels. As mentioned in the introduction, 
for typical Fermi moment u m 1 9 ' 2 1 Uf ~ 0.2 — O.Snm^ 1 , 
this would require a few nm precision in the height of 
the sample. It seems that, in the least, such sensitiv- 
ity to surface roughness would have to be eliminated in 
practice. 

In the Corbino geometry, the increase of the flux will 
transfer the charge from the inner to the outer surface, 
which can then measured. In the thermodynamic limit, 
this quantization is robust to the presence non-chiral 
states. Such measurement of a xy has been performed in 
2DEGs (Ref^) where both integer and fractional quanti- 
zations have been detected, and should be feasible in the 
3D topological insulators in the quantum Hall regime. 



1. Top horizontal surface (— <j>o < cj> < <j)o) 

For the top horizontal surface, r(cj)) — a/ cos cj) and 
— 4>o < <j> < cf>o where tan0o = b/a. For <f> in this range 
the metric for the surface is 



ds 2 = dx 2 + dx\ = dx 2 
= dx 2 + (da tan cj))' 



\4>) 



(dr{4>)\ 
\ # J 



dcf) 2 



(Al) 



In the Landau gauge of choice here, the two components 
of the vector potential are 



Ai = Bat&ncf), A 2 = 0. 



(A2) 



The physical spin is related to the Dirac Pauli matrices 
byl 



s = n x a 



(A3) 



where n is the normal to the surface. Using n = z and 
S1S2 = is3, we obtain 

si = -<r 2 , 82 = 01, S3 = 0-3- (A4) 
Thus, for —0o < 4> < 4>o the Hamiltonian is 



H = hvf 



,19 a 



1 d 



CT 2 - 



i da tan 4> 



h z a 3 . 
(A5) 



Ib = y j§ ■ This is clearly separable, and the eigenfunc- 
tions are planewaves in the x— direction. Its wavevector 
is set to be k. Letting p — y/2 (f^ tan cj> — ktB^j leads to 
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Appendix A: Derivation of the Hamiltonian, 
eigenstates and the matching conditions 

In this appendix we present the detailed steps which 
lead to the Eq.© as well as to the matching conditions 
which lead to the equation for the energy spectrum. We 
use the parametrization shown in Figll] 



h z -E 

(l hp 



The generic solution to these two coupled first order dif- 
ferential equations is 




0A6) 



aiD l ,(p) 
PiD v+1 [p) 



a 2 B v (-p) 
fcD v+l (-p) 



where D v (p) is the parabolic cylinder function^ 



D u (p) = 2»' 2 e- p ~ /i 

p n-h] 



r[(i - iy)/2i 1 1 V 2' 2' 2 



1-1/ 3 p 2 



y/2T[-v/2] 

and i-Fi is the confluent hypergeometric function 

t? 1 j. ^ 1 , a P , a ( Q + 1 ) P 2 , 
1 F 1 (a;6;p) = l + -- + - ( - TTy - + .... 



(A8) 



(A9) 



Unless v is a non-negative integer, the two solutions in 
the Eq. (|A7|) are linearly independent and D v {p) diverges 
as p — > — oo. These functions satisfy the relations 

§- p +\p) D »(p) = vD»-x{p), ( Al °) 

§- p -\p) D »(p) = -Du+i{p). (All) 

Since we are interested in taking the limit 6 3> (b, 
the solution near the outer surface must satisfy vanishing 
boundary condition as it moves closer to the 2D "bulk". 
This means that near the outer surface c\ l = and 



u k {4>)\ = ( (^ + Vz)Di (€ 2_ n 2 ) _ 1 {y/2{K-atBXi4>)) 
v k {4>) J ClH \ v^I>i(e»-^)(v^(K-atan^)) 



where the dimensionless lengths and energy scales are 



a=fz^=TZ^ = U B , e = -j-fj^, Vz = and 



w 9 hv F /e B ■ 

Similarly, near the inner surface Cx t u = and the so- 
lution has the form 

uk{4>) 

Vk(<t>) 

( -(e + ( e 2_^)_i (\/2(atan^- «)) V 

Cl 'M ^1(^2) (n/2(« tan 



2. Outer vertical surface (<j>o < 4> < t — 0o) 

r(</>) = 6/sin</> 

(is 2 = cLt 2 + e?a; 2 = <ix 2 + (<i6 cot </> 
=>■ x\ — x, x 2 = a + b — 6 cot <j). 

Ai = Bb, A 2 = 0. 
The physical spin operators (up to H/2) are 
Sl = c 3 , «2 = cr 2 , S 3 = -<7i. 
For </>o < < 7r — 4>o the Hamiltonian is 



(A17) 



H = Hvf 



Id eB\ 1 d 

fl|Tn ~ ° - °"2- 



96 cot i 



h z o\. 
(A18) 



and the eigenfunctions are 

f ] _ c V e2 -( K -^-^) 2 ' 9cot ^ i 
V M<t>) J 2,1 V 



e 



C2.2e 



_i v / e 2_ (K _ / 3„^ ) 2 /3c 



(W9) 



where 



a± 



K -/3-77^Ve 2 -( K -/3- % ) 2 . (A20) 



3. Bottom horizontal surface (7r — 0o < < tt + 0o) 

r((f>) — — a/cos</> 

ds 2 = dx\ + dx\ = dx 2 + (datan</>) 2 (A21) 
x a = x, x 2 = 2(a + b) + cit&ncj). (A22) 

Ai = -Batan0, A 2 = 0. (A23) 

Sl = 02, s 2 = -ffi, s 3 = cr 3 . (A24) 
^jl^r 7T — </io<0<7r + 0o the Hamiltonian is 

H = Hvf 



18 . o 

B 



CTll 7fe + F tan0l+ " 2 - 



i a 



i da tan 

(A25) 



where we included a different electrical potential on the 
bottom surface V g . 

Near the outer vertical surface the solution on the bot- 
tom horizontal surface is 



Uk(4>) 

Vk{(t>) 



(A26) 



( V2D i ((e _„ s) 2_ J7i) (V2 (k + a tan 0)) 

C3 ' fl \ (e - % - f 9 )£»| ((e _ I/s) 2_ j; 2 ) _ 1 (a/2 (re + a tan (^)) 

where i/ g = V g / {Tivf / ^b)- 

Near the inner vertical surface the solution on the bot- 
tom horizontal surface is 



(A14) 


( ~ 


(A15) 


1, Vk(4>) , 


(A16) 


C3,L , 



(A27) 



V / 2 J Di(( £ _„ s) 2_ J ,2 ) (-V^(K + atan0)) 
e + »y a )Di ((e _ 1/s )2_,2 ) _ 1 (-- s/2(« + atan^ 

4. Inner vertical surface (-7T + </>o < (f> < 2n — <j>o) 

r(4>) = —bj sinijf) 

ds 2 = dx\ + dx\ = dx 2 + (db cot (j>) 2 (A28) 
=S>xi = x, x 2 = 3(a + 6) - 6cot0. (A29) 



Aj = -56, A 2 = 0. 

Sl = —03, S 2 = CT 2 , S 3 = Oi. 

For 7r + 0o < < 27r — 4>q the Hamiltonian is 
15 b \ Id 



H = hvf 



i 96 cot 



(A30) 
(A31) 

/l z (Tl. 

(A32) 



The eigenfunctions are 

V "fcW J ' v e 

+ C 4 i2 e-V^ 2 -(«+^) 2 ^-^ C a - (V33) 



where 



a± = k + P + r] z ± i^e 2 -{k + /3 + r) z ) 2 . (A34) 



5. Matching conditions 

As discussed in the main text, we require the continuity 
of the wavefunctions near the outer surface where k w p. 
Therefore, we must have 



Uk{Tl ~ (j) { +) ) 



(A35) 



^-<f|)V (A36) 



where (±)o = O ± 0+. 

Using the wavefunctions determined in the Appendix, 
the above set of four linear equations in four unknowns 
translates into 



e ec 2 ,i 



e ec 2 ,2 



S+Cl,iJ 

A,+C3,fl = e ie - a ec 2 ,i+e~ ie - a ec 2 , 2 



c 2A +e- l0 - a (uj- -i8-)c 2 ,2 



i0-)c 2 ,2 
(A37) 



where 
At = 

A± = 

w± = 
0± = 



±(e + Vz )D He2 _„,)_! (V2(-/? ± «)) (A38) 

^i (e2 _ r) ,)(V2(-/?±«)) (A39) 

±(e - i/ fl - ((e _ I/a) 2_ r? 2 ) _ 1 (\/2(-/3 ± «)) 

(A40) 

^|(( £ -, 9 ) 2 -^) (^2(-/3 ± «)) (A41) 
k ± /3 ± 7? 2 (A42) 

^/e 2 (A43) 



Near the inner surface, where k w — /?, we must have 



V Wfe(7T + 4" ) ) 

/ u fe (2^-^ +) ) 
^(2tt-^ +) ) 
This translates to 



u k (ir + (f) ( Q +) ) 
v k (ir + (f) ( Q +) ) 



(AAA) 



(A45) 



Bg-C?,,L 

A-Ci, L 
B-d, L 



e lS + a (u+ + i6 + )c 4A + e- w + a {w+ - i6 + ) Ci , 2 



e- lf> + a (Lu+ + i0+)c4,i + e ie + a (u; + - i0+)c4,: 



ec 4 .i + e iW + a ec 4 , 2 



(A46) 



1 D.-H. Lee, Phys. Rev. Lett. 103, 196804 (2009). 

2 V. T. Dolgopolov, A. A. Shashkin, N. B. Zhitenev, S. I. 
Dorozhkin, and K. von Klitzing, Phys. Rev. B 46, 12560 
(1992). 

3 L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 
106803 (2007). 

4 L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 

5 D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, 
and M. Z. Hasan, Nature 452, 970 (2008). 

6 D. Hsieh, Y. Xia, L. Wray, D. Qian, A. Pal, J. H. Dil, 
J. Osterwalder, F. Meier, G. Bihlmayer, C. L. Kane, et al., 
Science 323, 919 (2009). 

7 Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. Ban- 
sil, D. Grauer, Y. S. Hor, R. J. Cava, et al., Nature Physics 
5, 398 (2009). 

8 Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S.-K. Mo, 
X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, et al., 
Science 325, 178 (2009). 

9 X. L. Qi and S. Zhang, Physics Today 63, 33 (2010). 

10 J. Moore, Nature 464, 194 (2010). 

11 M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 
(2010). 



12 R.-L. Chu, J. Shi, and S.-Q. Shen, Phys. Rev. B 84, 085312 
(2011). 

13 Z. Yang and J. H. Han, Phys. Rev. B 83, 045415 (2011). 

14 Y.-Y. Zhang, X.-R. Wang, and X. C. Xie, ArXiv e-prints 
(2011), 1103.3761v2. 

15 K. Nomura, S. Ryu, M. Koshino, C. Mudry, and A. Fu- 
rusaki, Phys. Rev. Lett. 100, 246806 (2008). 

16 A. A. Taskin and Y. Ando, Phys. Rev. B 80, 085303 (2009). 

17 J. G. Checkelsky, Y. S. Hor, M.-H. Liu, D.-X. Qu, R. J. 
Cava, and N. P. Ong, Phys. Rev. Lett. 103, 246601 (2009). 

18 J. G. Analytis, R. D. McDonald, S. C. Riggs, J. H. Chu, 

G. S. Boebinger, and I. R. Fisher, Nature Physics 6, 960 
(2010). 

19 J. Xiong, A. C. Petersen, D. Qu, R. J. Cava, and N. P. 
Ong, ArXiv e-prints (2011), 1101.1315. 

20 D.-X. Qu, Y. S. Hor, R. J. Cava, and N. P. Ong, ArXiv 
e-prints (2011), 1108.4483. 

21 C. Briine, C. X. Liu, E. G. Novik, E. M. Hankiewicz, 

H. Buhmann, Y. L. Chen, X. L. Qi, Z. X. Shen, S. C. 
Zhang, and L. W. Molenkamp, Phys. Rev. Lett. 106, 
126803 (2011). 

22 A. Pnucli, Journal of Physics A: Mathematical and General 



27, 1345 (1994). 

M. Stone and P. M. Goldbart, Mathematics for Physics: A 
Guided Tour for Graduate Students (Cambridge University 
Press, 2009), 1st ed. 

J. H. Bardarson, J. Tworzydlo, P. W. Brouwer, and 
C. W. J. Beenakker, Phys. Rev. Lett. 99, 106801 (2007). 
P. M. Ostrovsky, I. V. Gornyi, and A. D. Mirlin, Phys. 



Rev. Lett. 105, 036803 (2010). 

M. Sitte, A. Rosch, E. Altman, and L. Fritz, ArXiv e-prints 
(2011), 1110.1363vl. 

E. Merzbacher, Quantum Mechanics (John Wiley & Sons, 
New York, NY, USA, 1998), 3rd ed. 



